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Abstract 

We employ the Polyakov world-line path-integral version of QCD to iden- 
tify and resum at leading perturbative order enhanced radiative gluon con- 
tributions to the Drell-Yan type {qq pair annihilation) cross-sections. We 
emphasize that this is the first time that world-line techniques are applied to 
cross-section calculations. 
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I. INTRODUCTION 



The path, rather than the functional, integral casting of a relativistic quantum system 
has a long history, going back to Fock [1], Feynman [2], and Schwinger [3]. In the course of 
time it has received substantial contributions from several authors (see, e.g., Refs. [4-9] to 
name just a few). However, it was not realized until recent developments in string theory 
in the context of effective actions (see [10-12], and [13] for a recent review) that first- 
quantization methodologies in High Energy Theory can compete with second-quantization 
ones. Of particular interest to us here is the Polyakov world-line path integral [14], which 
employs world-line paths weighted by a spin factor with the aim to describe the propagation 
of particle-like entities in Euclidean space-time. In fact, Polyakov's intention was to use 
this construction as a simple prototype for discussing string quantization. Hence, for his 
purposes it was sufficient to consider the simple case of a free, spin-1/2 particle-like entity. 
Motivated by this, two of the present authors [15] explored the possibility of transcribing the 
matter, spin-1/2, field sector of a gauge theory into a Polyakov world-line path-integral form. 
In these works, it was established, for both Abelian and non-Abelian gauge systems, that 
this is, indeed, possible with the spin-factor making explicitly its entrance in the resulting 
expression, while the dynamics enters through a Wilson line (loop) factor defined on each 
given path. 

Due to the Gaussian character (with respect to the Dirac fields) of the fermionic sector 
of physically relevant gauge-field theories, the aforementioned transcription into a Polyakov 
world-line path-integral refers to the full system. This means that one's way of thinking 
should be readjusted to the idea that the second-quantization formalism, associated with 
the field theoretical mode of description, can be replaced by a new, but equivalent, structure 
that is based on space-time path integrals. The quantities of central importance defined 
on these paths are then the spin-factor and the Wilson line (loop), the latter becoming an 
indigenous element of the theory, as it enters at the level of its definition. 

In a number of papers see, for instance, Refs. [6,15-17] - we have employed the path- 
integral casting of either QED or QCD, to study infra-red (IR) factorization and ensuing 
behavior of Green's functions and amplitudes in a resummed perturbative context at the 
two-, three- and four-point function level. Roughly speaking, the aforementioned isolation 
of the long-distance physics in these theories emerges through the ability to identify a spe- 
cial set of space-time paths having a very simple geometrical profile which is shared, in a 
restricted (but directly relevant to the physics of the process) neighborhood, by each and 
every contour entering the path integral. In a Euclidean space-time context, the single (mul- 
tiplicative) renormalization constant, carried by this special family of paths, automatically 
factorizes out [18] their contribution to amplitudes/cross-sections, given that it also accom- 
panies the rest of the paths. The more complex geometrical structure of the latter, simply 
implicates additional ultraviolet (UV) singularities which can be absorbed into conventional 
wave function and coupling constant renormalizations. This clean, geometrically based, ar- 
gument, which singularly underlines the world-line description, will be further elucidated 
through the main exposition in the sections to follow. Minkowski space subtleties, associ- 
ated with the light cone, which are encountered in the particular processes under study, will 
require separate attention. 

Perhaps the most important accomplishment of this paper is that it extends world-line 
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techniques to cross-section calculations for the first time. To be sure, the situation presently 
considered refers to the hypothetical situation where the scattering process involves quarks 
and associated gluon radiation without reference to hadrons. It does, nevertheless, fall 
within the spirit that marks our approach to IR issues in QCD [16,17]: Once off-mass shell 
IR protection is employed - by an amount that exceeds Aqcd - one actually tests how far one 
can go by remaining strictly within the confines of QCD before attempting to make contact 
with real hadrons. Granted the opposite route, from hadrons to quarks and gluons, via the 
use of quantities like structure/fragmentation functions, the employment of tools such as the 
operator product expansion, etc., constitutes a more realistic procedure for investigating the 
same physical problems. On the other hand, an effort which bases its considerations on a 
fundamental theoretical framework in order to arrive at "cross-sections" does present merits 
and interests of its own, an example of which will be presented below. In this context, the 
philosophy underlying our approach to the IR domain of QCD is closer in spirit to the one 
articulated by Ciafaloni in [19], the only difference being that we shall keep a more pragmatic 
(and less ambitious) course by focusing our attention on cross-section expressions. 

Letting these comments suffice for an introductory exposition, we now proceed to dis- 
play the organization of the paper, which is as follows. In the next section, we exhibit the 
world-line expression for the full fermionic Green's function and subsequently employ it to 
construct corresponding expressions for DY-type QCD amplitudes/cross-sections. Section 
III furnishes, with the aid of an Appendix, our basic calculations associated with one virtual 
gluon exchanges for the specified special set of trajectories. The resulting expression explic- 
itly reveals the threshold enhancement factor, whereas the task of virtual gluon resummation 
is performed, via the aid of the renormalization group, in Section IV. Section V deals with 
the resummation of contributions from real gluon emission. Finally, in the last section, we 
further discuss our results and present our conclusions. 



II. BASIC WORLDLINE EXPRESSIONS FOR AMPLITUDES AND 

CROSS-SECTIONS 



Consider the full two-point (fermionic) Green's function in the presence of an external 
gluonic field. The Polyakov path-integral expression, in Euclidean space-time, 
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displays the basic world-line features pertaining to this quantity. Here, and below, P denotes 
the usual path ordering of the integrals. The first thing to point out is that a given path 
of the matter field quantum, starting at x and ending at y between respective "proper- 
time" values and T, also enters a Wilson line factor. The latter, being the sole carrier 
of the dynamics, separates itself from the rest of the factors in the path integral which are 
associated with geometrical properties of paths traversed by spin-1/2 particle entities. The 
most notable such quantity is the so-called spin factor [14], Pexp Ui/A) dt a • uj , where 
oj^ u = (T/2)(x fJ ,x l/ — Xfj,x u ), accounting, in a geometrical way, for the spin-1/2 nature of 
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the propagating particle. Accordingly, our perturbative expansions should be perceived of 
in terms of (Euclidean) space-time paths involving a "proper time" parameter and not in 
terms of Feynman diagrams. As it turns out [20], in the perturbative context, the structure 
of matter particle contours, entering the path integral, is determined by the points, where a 
momentum change takes place, i.e., points where a gauge- field line (real or virtual) attaches 
itself on the (fermionic) matter-field path. The almost everywhere non-differentiability of 
these contours, is residing precisely at these points. A major effort, in this paper, will 
be devoted to the extension of the world-line formalism to expressions for cross-sections 
corresponding to the particular processes of qq annihilation. 

From the world-line point of view, the process we intend to study involves fermionic 
matter particle (quark) paths that commence at x and end at y, being forced to pass through 
an intermediate point z, where a momentum transfer Q takes place. This means that the 
Green's (vertex-type) function we shall be dealing with has the following form (r^ denotes 
some Clifford-Dirac algebra element) 
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where 
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It is especially important to realize that in our approach off-shellness is naturally parame- 
terized in terms of the finite size of the matter particle contours and realistically accounts 
for the fact that quarks reside inside a hadron (m can be viewed as an effective quark mass). 
Going over to momentum space, we write 
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where C z denotes a generic path forced to pass through point z, at which the momentum Q 
is imparted. 

For a process of the type q + q — > lepton pair + X the "amplitude" expression reads 
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with the second, comprehensive, expression to be understood having recourse to Eq. (4). 

For the cross-section, we need to employ the following quantity, which we implicitly 
display in Minkowski space-time after straightforward adjustments, 
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where P denotes anti-path ordering. Even though not explicitly displayed, the cross-section 
acquires a path-integral form, which has the following characteristics: 

1. Paths C z and C z ' are forced to pass through points z and z' , respectively, where the 
momentum transfer occurs (see Fig. 1). The distance b = \z — z'\ serves as a measure 
of how far apart the two conjugate contours can venture away from each other and 
will be referred to as the impact parameter. 

2. The traversal of C z ' is made in the opposite sense relative to C z . If now, we let the 
two paths join at one end by using translational invariance, while we allow the other 
two ends of the contour to close at infinity, then we obtain the formation of a Wilson 
loop. 

3. Under these circumstances, the Wilson loop formation guarantees the gauge invariance 
of the expression for the cross section. 

On the other hand, by keeping the contour lengths finite, but very large, thereby placing 
the quarks off- mass-shell, gauge invariance will still continue to hold to the order of approxi- 
mation we employ in our computations, given that the off-mass-shellness serves at the same 
time as an IR cutoff. 

Up to this point our considerations have been centered around the geometrical profile 
of the paths entering the world-line casting of QCD, the main conclusion being that, for 
the process considered, the relevant contours entering the path integral are marked by a 
characteristic point, where a momentum transfer is imparted and that they are open for the 
amplitude and closed (or almost so) for the cross-section. Armed with this information, we 
now turn our attention to the Wilson factor which contains all the dynamics of the given 
process. The obvious task in front of us is to assess its implications once the gauge fields 
are quantized, i.e., once the Wilson factor is inserted into a functional integral weighted by 
the exponential of the Yang-Mills action. We display the quantity of interest as follows 
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In the above expression, {• • -}a signifies the expectation value with respect to the gauge 
field functional integral which, in this work, will be considered in the context of perturbation 
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FIG. 1. Illustration of two conjugate contours C and C entering the world-line path integral, "talking" 
to each other at points z and z', where the momentum transfer for the physical process takes place. The 
distance \z — z'\ is referred to as the impact parameter. 

theory. Note in the same context that a virtual gluon attaching itself with both ends to the 
fermionic worldline, entering the amplitude, corresponds to a correlator between a pair of 
gauge fields originating from the expansion of the Wilson factor. On the other hand, for an 
emitted "real" gluon from the fermionic line, the correlator is between an "external" and 
a Wilson-line gauge field. 1 The overall situation is depicted in Fig. 2. At the cross-section 
level, now, "real" gluons are integrated with respect to "propagators" linking together the 
two conjugate contours, while their polarization vectors are summed over (cut propagators). 
This is precisely what (• • •) signifies in the last equation, as it brackets both Wilson line 
factors. 

This marks a crucial difference to conventional approaches (for example, [21-23]), wherein 
the Drell-Yan process is discussed in a context where IR factorization is based on the eikonal 
approximation for soft amplitudes. Wilson loop expectation values, entering this scheme, are 
evaluated along contours corresponding to classical trajectories - along with a segment which 
lies on the light cone, introduced in order to secure gauge invariance. In our case, by contrast, 
Wilson contours are built in at a foundational level, being themselves an integral part of 
the description of the full QCD. Accordingly, factorization properties for us are integrally 
connected with the renormalization properties of Wilson loops studied in the more general 
context of Refs. [18]. In the light of the above remarks, let us proceed to display the first- 
order (in perturbation theory) expression for W, which receives contributions from virtual 



One will, of course, also encounter correlators that involve gauge fields from the non-linear terms 
of the Yang-Mills action. These, however, do not enter the leading logarithmic considerations 
relevant for our considerations. 
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gluons, viz., those attached at both ends of either the world-line contour C z or C z ' , as well 
as from "real" gluons linking these contours to each other (cf. Fig. 2). This expression reads 

W (2) = Tr / - g 2 C F £ dt, £ dt 2 9 (t 2 - x» (t 2 ) x v (h) (x(t 2 ) - xfa)) 

- g 2 C F [ T dt[ [ T dt' 2 9 (t[ - Q & (t' 2 ) x v (t[) D, u (x(t' 2 ) - x(t[)) 

Jo Jo 

- g 2 C F £ dt £ dt' x(t) ■ x (0 D cut (x(t) - x(t)) + O (g A ) . (8) 

It becomes obvious from their structure that the first two non-trivial terms correspond to 
virtual gluon contributions - one per conjugate branch while the third one is associated 
with "real" gluon emission. Finally, concerning the gluon propagators entering the above 
equation, we shall be employing their Feynman-gauge form without loss of generality due 
to gauge invariance. In particular we have, in D-dimensions, 

4-d f d ° k e ~ ik ' x 1 ( 2 \( 4 ~ D)/2 r(D/2-i) 

whereas 

D cut (x)=^- D J -^2vr5(g 2 )^(g°)e-^ 



From here on and for the sake of notational simplicity, we shall simply write D{x) instead 
of D cnt (x). 

As already established by other methods, the perturbative expansion (8) is plagued by 
large threshold logarithms leading to the need for factorization and resummation. This is 
precisely the task we are about to undertake within our framework. 

III. FIRST-ORDER VIRTUAL GLUON CORRECTIONS IN THE VICINITIES 

OF POINTS z AND z' 

The space-time mode of description of the Polyakov world-line formalism puts us into 
the position to promote the following argument: The point z (or z'), where the momentum 
transfer Q is imparted, marks the presence of a neighborhood around it, no matter how 
infinitesimal in size this might be, whose geometrical structure is shared by all fermionic 
paths entering the path-integral. Specifically, there will be a derailment (cusp formation), 
whose opening angle will be fixed unambiguously, since it is determined by the momentum 
transfer. It follows that the contributions to the amplitude and cross-section from the imme- 
diate vicinities of each of the two cusps is a common feature of all contours and eventually 
factorizes. In this section we shall determine the first-order perturbative term corresponding 
exactly to this factor. 

Consider now the neighborhood of point z on the contour C z . Expanding around this 
point, we write 
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x(s') = z 



z = x(s) 



FIG. 2. Virtual gluon radiative corrections of various sorts and "real" gluon lines with their ends 
attached on each of the two depicted contours at the cross-section level. 



x »(t) = x"{s) + (t- s)^(s ± 0) + . . . 



(11) 



with v M = x^(s — 0) and v ' M = x^(s + 0) being entrance and exit four- velocities, respectively, 
with respect to z. 

Adjusting our notation by re-parameterizing the contour so that the zero value is assigned 
to point z, the relevant quantity to compute, to first perturbative order, becomes 
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It is clear that the above expression corresponds to the first term in Eq. (8), which monitors 
a virtual gluon exchange occurring on contour C. 

From the above considerations it follows that the main contribution to each double in- 
tegral comes from the common limit t±, t 2 — > 0. Suppose now, the other limit is to be 
determined by demanding that its contribution to the integrals is of vanishing importance. 
Then, such a requirement automatically isolates those contours, whose only significant ge- 
ometrical characteristic is that the four-velocities to approach and depart from point z are 
fixed, denoted by v M and v'^ respectively. The same, of course, happens for point z' , but 
in the reverse order. This justifies the subscript S in Uc\s, which stands for "smooth". 
Let us also observe that the omitted terms in Eq. (12) will contain negative powers of a, 
whose dimension in the denominator is (mass) 2 . Neglecting their presence means that o 
should be very large in magnitude and hence it should be related to an IR cutoff, i.e., 
a ~ A" 2 , A > Aqcd- 

In Euclidean space-time, now, every path will share the geometrical structure we are 
focusing on in some neighborhood of the point z or (z'), no matter how close to these points 
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one has to come. At the same time, the UV singularities, exhibited by this restricted set of 
paths, will entail expressions that solely depend on the two four-velocities and the opening 
angle. Paths of more complex geometrical structure, on the other hand, will certainly 
exhibit these UV singularities plus additional ones. 2 It follows that - in Euclidean space- 
time at least - the restricted set of trajectories, by exclusively carrying the corresponding 
(multiplicative) renormalization constant, factorizes from the rest of the expression for the 
amplitude and/or cross-section. In a Minkowski space-time context, which will be considered 
next, we should anticipate the existence of additional contributions to Eq. (12), due to the 
light-cone structure that cannot be assigned to each and every contour and, therefore, cannot 
be factorized. Let us, then, go over to Minkowski space-time, where we have two distinct 
possibilities for defining an infinitesimally small neighborhood around z. The first one, to 
be labelled (a), reads 

(x-x'f = C(e 2 ), with ~ v'^ for all \i, (13) 

where e(< Q~ v ) is a small length scale. The second alternative, to be labelled (b), can be 
typically represented by 

{x - x') 2 =0{e 2 ) with \v - v'\ 2 = 0{\ 2 ) but (y+ - v' + ) ~ 0(Q) and (u_ - v'_) ~ O 

^(v + -v' + )(v--v'_) = 0(\ 2 ) (14) 

that is equivalently effected via the condition v + ^> v' + , i>_ ~ v'_. All in all, there are four 
different configurations: + <-> — and prime^no-prime entering this case. 

We denote case (a) as "uniformly soft" , given that the considered gluon exchanges take 
place in a neighborhood whose smallness pertains to all directions. Case (6), on the other 
hand, will be referred to as "jet" since gluon emission occurs under circumstances, where 
entrance and exit four- velocities differ from each other significantly along one or the other of 
the light-cone directions. Particular implications stemming from this, purely Minkowskian, 
case as far as the factorization issue is concerned, will be considered later on. 

Let us commence our calculations by taking up the first 0(g 2 ) term entering the right 
hand side of Eq. (12). Since this only involves the branch of the contour C z entering point 
z, we obtain the same expression regardless of whether or not a uniformly soft or a jet 
configuration is being considered. It reads 

h = dh dt 2 6 (t 2 - ti) vWD^ (vt 2 - vh) 

J —a J — a 

= -si(-Vi;) r (T- 1 ) -32- g/ 2 ' (15) 

where 3 L x = cr\v\. The second term has the same structure as the first one (it involves the 
exiting branch of C z ) and therefore produces a similar result: 



2 Actually, the standard UV singularities of perturbative field theories associated with (3- 
functions, coupling-constant and wave-function renormalization, pertain to almost everywhere non- 
differentiable paths. 

3 Note that v has dimensions of mass as our "time" parameter a has dimensions of (mass) -2 . 
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N (4-D)/2 ^ ( D 



'' = -wt*M) r {l- l )— 3 2^DT2 (16) 

with L 2 = cr\v'\. 

A couple of remarks are in order at this point. First, even though the length scales L\ 
and L 2 are both large, being proportional to a, they will be of the same order of magnitude 
for case (a), whereas for case (b), one scale will be negligible in comparison with the other. 
Accordingly, the total expression to the amplitude for the uniformly soft contribution will 
be twice as large as that of the jet-like one. This being said, we shall denote the dominant 
length scale by L(~ L\ and/or L 2 ), when it enters our final expressions, and set it equal 
to 1/A, recognizing that it is of the same order as the IR cutoff. Second, in order to avoid 
the double counting resulting from the fact that each branch has been "cut-off" at distance 
L away from z, where gluon emission occurring at the endpoints will be offset by a similar 
one, but opposite in sign, from that portion of the contour that continues to stretch out 
to infinity, the final expressions for the end-point singularities should be multiplied by a 
factor of 1/2. Equivalently, one might think of this compensation as actually identifying 
the missing energy of the gluon emission at the extremities of the path with the off-mass- 
shellness. In fact, this is what we have been implying all along when claiming that finite 
contours signify off-mass-shellness. 

Turning our attention to the contribution resulting from a virtual gluon exchange from 
the entrance to the exit branch, with respect to z, we consider the quantity 



xv ■ v' j o dh j o dt 2 (tlv 2 + t 2 2 v' 2 + 2ht 2 v ■ v' - i0 + ) . (17) 

For case (a) it assumes the form (recall that v ■ v ' is negative) 



h 



(a) 




1-0D/2) 



(18) 



As shown in the Appendix, one then determines (7e is Euler's constant) 

^ = ^coth 7 ^ + ^ 7 coth 7 ln (^e 2 +^ , (19) 
where cosh 7 = w = — , v ;f., > 1 . 

' \v\\v'\ — 

In all of the above expressions, as well as in those that will follow, we have ignored: 
(i) all imaginary terms that will drop out when contributions (for virtual gluons) from the 
conjugate contour are taken into account and (ii) finite, //-independent terms that will cancel 
out when real gluon contributions to the cross-section are included. 

Collecting all terms, we deduce, for the "uniformly smooth" part, 

I? + 4 a) + # = ^(7Coth 7 - + ^(7Coth 7 - 1) In (^^A ■ (20) 
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Concerning the "jet" part of the computation, we only need to consider 1^ because 4 
the expression for l[ b ^ + 1^ is simply one half of that of l[ a ^ + A typical term entering 
4 b) (u_ > v' + ) is 

1 / o\ (4— -D)/2 / D \ , t" „ r „ /.2 2 , „. , / .„ \M.D/ 2 ) 



/(») = JL ( 4 -^)/ 2 r - 1) w jT dti jT dt 2 (*y + 2t 1 t 2V ■ w ' - *o + ) 1 



(21) 



whose computation suffices to furnish each of the other three terms as well. 
It is shown in the latter part of the Appendix that one obtains 

lf = is? (^H + ^ ln2 + const - < 22 > 

It is duly observed that the singularity structure of the above expression is 7-independent. 
In fact, the "jet" configuration is a Minkowski-space feature and is connected to "gluon" 
emission in the + or the - light-cone coordinates direction. This result is in accord with 
Wilson loop studies in Minkowski space, wherein the relevant contour lies partly on the light 
cone [22]. Subtracting the pole terms in the MS scheme, we arrive at the finite part of the 
overall result. For the uniformly soft contribution, in particular, we get 

(l{ a) + 4 a) + 4 a) ) fin = JL ( 7 coth 7 - 1) In (^j , (23) 

while the jet contribution reads 



167T 2 A 2 

where we have set A 2 = 4A 2 e~ 2lE . The above relation takes into account all four different 
configurations contributing to I%\ 

Gathering all terms, we arrive at the following overall result for the second-order contri- 
bution stemming from contour C z 



c ' s 2tt 



( 7 coth 7 -l)ln(^]-^ln(^]+ln 





(25) 



A similar result is obtained also for contour C z ' . 

Noting that 7cotli7 = In (Q 2 /m 2 ) (for Q 2 ^> m 2 ), with Q 2 = {p + p') 2 , we recognize 
that the well-known perturbative enhancements occurring as Q 2 — > 00 are associated with 
the eikonal-type trajectories upon which our present calculations have been based. One, 
now, realizes that these trajectories define threshold conditions, with respect to the given 
momentum exchange Q, for the process under consideration, since they leave no room for 
space-time contour fluctuations. In the following section, we shall treat the resummation of 
these enhanced contributions to leading logarithmic order. We shall, furthermore, identify 
a correction factor associated with those terms in Eq. (25) not involving the enhancement 
factor ln(Q 2 /m 2 ). 



4 Recall the remark following Eq. (16). 
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IV. RESUMMATION OF ENHANCED CONTRIBUTIONS FROM VIRTUAL 

GLUONS 



The family of world-line paths to which the considerations in the previous section refer 
was used in order to deal with all (virtual) single-gluon exchanges, consistent with the simple 
geometrical configuration of two constant four-velocities making a fixed angle 7 between 
them (in Euclidean formulation). Among these gluons there will be "hard" ones (upper limit 
Q) and "soft" ones (lower limit set by A). What is debited to the former and what to the 
latter group of gluons is, of course, relative. It is precisely the role of the renormalization scale 
fi, entering through the need to face UV divergences arising even for the restricted family 
of paths, to provide the dividing line. The corresponding renormalizat ion-group equation 
reflects the fact that the scale /i is arbitrary and that physical results do not depend on it. A 
straightforward application of this fact will enable us to resum the enhanced, virtual gluon 
contribution to the amplitude in leading logarithmic order, as well as to obtain a bona-fide 
correction term. 

To bring the above discussion into a concrete form, let us first consider a separation, 
good to order 1/Q 2 , of the cusp contribution, which can be factorized from the amplitude 
Uc, entering Eq. (7), on the basis of what has been determined so far. Hence, we write 



Uc=u ^(^yj^ +0 (^ (26) 



r(2) ■ a « 



with 

t^cusp = 1 - gC F (7 coth 7 - 1) In (j^J , (27) 

where we have normalized t/c iCUSp to unity for 7^0. The designation "cusp", above, refers 
to that factor of the soft sector, which recognizes the angle 7. The factor U c , on the other 
hand, includes both: (i) soft contributions - related to the dependence on the quantity fi 2 / A 2 
- and (ii) hard ones - depending on the quantity Q 2 / fj, 2 . 

It is convenient to take the logarithmic derivative of Eq. (26) with respect to Q 2 : 

^ 1,1 Uc = 5W 1,1 u »~ + l^ hl0c + o (£) ' (28) 

The /^-independence of physical results leads to the renormalization group equation whose 
ultimate justification has to do with the multiplicative renormalization of the soft (cusp- 
angle dependent) factor. Indeed, the latter is detached from collinear emission and totally 
complies with the Euclidean space-time properties of Wilson loops for which the results of 
Refs. [18] fully apply. Specifically, we write 

\nU c = - f d rn£/ c ,cusp = r cusp (a s ) (29) 



d In fj, d In Q 2 d In // d In Q 2 

with T cusp to be read off from Eqs. (20)-(22) and (12) 

. 2 



T cusp (a s ) = ^C F + 0(a 2 ). (30) 

71 
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From the second leg of Eq. (29), one obtains 

d 



d\nQ'- 



■InU, 



C,cusp 



x 2 2t cusp 



(31) 



which, in turn, gives 



d 



d\nQ 2 
where we have defined 

ir 

2 



In 



rQ 2 dt. 



Uc = - I ^r cusp [a s {t)\ + T a s {Q 2 ) 



a s (Q 2 ) 



d 



d\nQ 2 

Combining the last three equations, we have 



(32) 



(33) 



m 2 =Q 2 



d 



\nU, 



dt. 



c 



d\nQ 2 

where, in terms of InUc, we write 

1 



- J ^ r ™s P Mt)} + -r [a s (Q ' 



-r 

2 



a s (Q< 



d 



lnU c \ 



\ 2 =Q 2 



dlnQ 2 

Setting /i 2 = Q 2 in Eq. (26), we are led to the identification 
1 



(34) 



(35) 



-r 

2 



a s (Q 2 ) 



dlnQ 2 
To second order we have 



d . T . (Q 2 Q 2 \ , d - / Q 2 ^ 

In f/c,cusp ( , ) I a 2 =Q 2 + "TT^ lnC/ c I 1, — 



F (2) 



« S (Q 2 ) 



3Gp 
2~T 



dlnQ 2 



a s (g 2 )+(9(a 2 ). 



A 2 



(36) 



A2=Q2 



(37) 



Gathering our findings, we obtain our final, resummed result corresponding to the contour 
C. It reads 



J f Q2 dt 

Uc=exp rU 2t 



Q 2 

In— r cusp (a s (t)) - T(a s (t)) 



Uc,oMQ 2 )). 



(38) 



One notes that the second ("correction") term in the square brackets is associated with 
collinear emission (cf. Eq. (36)). Finally, the factor U C fi{ct s {Q 2 )) represents input from 
initial conditions at the QCD level. Clearly, the conjugate-contour term U\C Z ') can be 
treated in a completely analogous fashion. 
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x(s') = z' 



z = x(s) 



FIG. 3. Neighborhoods of respective points on two conjugate contours, where the momentum transfer 
takes place, and associated four-velocities. 



V. RESUMMATION OF ENHANCED CONTRIBUTIONS FROM REAL GLUON 

EMISSION 

We shall now turn our attention to real gluons and attempt to factorize cross-section 
contributions from neighborhoods around points z and z' . Note that this time we have to deal 
with gluons which connect two "opposite" neighborhoods while crossing the unitarity line 
(this situation is depicted in Fig. 3). The relevant scale promptly entering our considerations 
is the impact parameter b = z — z', which must be eventually integrated over in order to 
get the physically measurable cross-section. Naturally, the short-distance cutoff in this 
integration will be provided by the (length) scale 

For the eikonal-type family of paths, and in first-order perturbation theory, the relevant 
quantity on which our quantitative considerations are to be based, i.e., the counterpart of 
Eq. (12), is given by 



U%L = l + g 2 C F 



cc,s 



/U rV 
dti / dt 2 v ■ v D (t±v — t 2 v + b) 
-a J —a 

+ r dti r dt 2 v> ■ v' D (tW - t 2 v' + b) 
Jo Jo 

/0 rcr 
dti / dt 2 v ■ v' D (tiv — t 2 v' + b) 
-a Jo 

+ J dti J dt 2 v' ■ v D (tii/ - t 2 v + b) , 



(39) 



where the bar denotes four- velocities for the conjugate contour and the subscript cut is 
henceforth omitted. 

(2) 

To identify the leading behavior of (7^ s , with respect to b, we shall consider first 
the situation corresponding to b = 0. The subsequent emergence of UV divergences, once 
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handled through dimensional regularization, will introduce a mass scale // that will be 
bounded from below by an IR cutoff A and from above by the (mass) scale 1/6. The resulting 
renormalization group equation will facilitate the resummation of the leading terms, just as 
in the virtual-gluon case. 

Let us start with our quantitative considerations by looking at the term 

Ji(6) = v ■ v / dh / dt 2 D {t x v - t 2 v + 6) (40) 

J— (j J —a 

with v 2 = v 2 = —v ■ v (see Fig. 3). 

Setting 6 = and using the expression for the cut propagator as given by Eq. (10), we 
obtain 

m - -± (-^?) <"* F (f - l) ^ ^ [1 - (2- - „] ,4!) 

which actually coincides with what one would obtain if the regular propagator was substi- 
tuted. The significance of this occurrence is that it leads to the same anomalous dimensions 
for the running of the real-gluon contribution to the cross-section as for the virtual part. 
This fact can be immediately verified via a direct comparison with Eq. (15). 
Isolating the finite part of the above expression, we write 

=-<>(£)■ (42) 

It is trivial to see that the same result holds also for J^'*" 1 . 
We next turn our attention to the term 

J 3 (6) =vv'[° dh f 7 dt 2 D (hv - t 2 v + 6) . (43) 

J -a JO 

Its computation will concurrently allow us to determine J±(b), which corresponds to the 
exchange primes no-prime in the expression above. 

Dimensionally regularizing the cut propagator, we then obtain 

j 3 (0) = ^(-^(-^r (f - 1) vv>l°J tl 

x £dt 2 (t\v 2 + t\v' 2 + 2v ■ v't x t 2 - i0+) 1D/2 . (44) 

Once again we record, by referring to Eq. (17), coincidence of the singularities and, by 
extension, of associated anomalous dimensions between virtual and real gluon expressions 
that contribute to the cross section. 

For the "uniformly soft" configuration the corresponding result is 



An 2 A 2 \2 J \v\\v 

i rl ( V • v' 



1-13/2 



x dtj dt 2 [t( + t 2 2 + 2t 1 t 2 ——-i0 + \ . (45) 
Jo Jo V \v\ \v \ I 
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Taking into consideration that 



\v\\v'\ 
12 



\v'\\v\ 



cosh 7 > we obtain 



jt\o) = 4 a \o) = ^(-^) (4 " D)/2 r (f - 1) cosh 7 

X Jo dtl I dt<2 ( t i+ t 2 + 2t i t 2cosh7-i0+) 



l-D/2 



whose finite part reads 



jW> fin (0) = jr- nn (0) = — 7 coth 7 In 



r (o),fin 



A 2 



/2" 



(46) 



(47) 



Turning now our attention to the "jet" configuration, we can actually go directly to 
J^ b \0), since j{ b \0) + .4^(0) furnishes half of the contribution of its uniformly soft coun- 
terpart, the reason being the same as the one given in the virtual gluon case. We thus 
have 



j(b) 
•J a 



(0) 



1 



4tt 2 



(4-D)/2 



D 



- 1 



V ■ V 



X 



Jo Jo 



v ■ v' 



l-D/2 



dt 2 t( + 1t x t 2 —— - iO 



(48) 



with an analogous expression holding also for j[ b \o). 

For the finite parts of the "jet" contribution, one obtains 



J 1 (6) ' fin (0) + J 2 (6),fin (0)+4 



,(6),fin 
J 3 



(0) + Jf ' fin (0) 



4tt 2 



In 2 





(49) 



Collecting our findings from the real-gluon analysis to the second-order level, we write 
for the finite contribution to the cross-section 



U. 



(2) 



cc,s 



7T 



12 



/ a' 

(7coth 7 - l)ln ( — 




/2 X 



- 2 



(50) 



At the same time, the singularity structure of the full expression for the cross-section 
entails a multiplicative renormalization factor, which is common to all "Wilson loop" con- 
figurations entering its description, but which is the only one that pertains to the family of 
eikonal-type paths under consideration. The reasoning is, of course, identical to the one given 
for the virtual gluon case. Therefore, the corresponding contribution to the cross-section 
factorizes and the same resummation procedure can be employed as for the virtual-gluon 
case. As already observed, the anomalous dimension is in both cases the same. There are, 
however, the following notable differences. First, the upper limit for the momentum of real- 
gluon emission is 1/b 2 instead of Q 2 . Second, there is a difference of sign, which becomes 
evident by comparing Eq. (25) with Eq. (50). Finally, no hard real-gluon emission enters our 
considerations - by definition. In this light, it is practically self-evident that the resummed 
expression for real-gluon emission becomes 



U, 



cc 



exp 




Q 2 

In — r cusp (a s (t)) - T (a s (t)) 



u, 



cc,o> 



(51) 



16 



where c = 4e _27E corresponds to the canonical choice [24]. 

We can now bring together real and virtual gluon results by referring to our generic 
expression for the cross-section as given by Eq. (7). First, we write 



W = (Tt {U\C Z ')U(C Z ))) = Uc, cusp Uc, cusp U c c, msp W + o (i 



(52) 



where the factor W contains both hard and residual soft contributions. 

Then, at the cross-section level, our threshold resummation of the virtual gluons reads 

rQ 2 dt 



UnU, 



C U C = cxp 



p t 



i Q2 r 
In — 1 

t 



cusp 



[a s (t)) - T(a s (t)) 



Thus, combining the above expressions, the final result reads 



W = exp I - / 

J c 



Q 2 dt 

c/b 2 t 



Q 2 

In— r cusp (a s (t)) - T(a s {t)) 



Wo, 



(53) 



(54) 



with r cusp and T given by Eqs. (30) and (36), respectively, an expression obtained before in 
[21], employing Wilson lines as a quantity attached to quark current operators. 



VI. CONCLUDING REMARKS 

In this paper we have applied first- quantization techniques to study threshold resum- 
mation of soft gluon radiation for DY-type of processes in QCD. We have addressed our 
efforts in an energy regime whose lower cutoff is high enough to justify an analysis in which 
reference to "gluons", as dynamical degrees of freedom, continues to make sense. 5 

First-quantization approaches to the study of relativistic quantum systems, involving 
either strings or world-line agents as their basic tools, have been employed by a number of 
authors during the last decade (the interested reader is recommended to the recent reviews 
[9,13]) as a viable alternative to the traditional second-quantization procedure, associated 
with their field theoretical casting. The main feature of this type of approach at the pertur- 
bative level, is that it allows a space-time description of the field theoretical system leading 
to expressions that accommodate a host of Feynman diagrams at once. On the other hand, 
it can play a crucial role in the study of non-perturbative effects [9]. 

The chief issues of the present methodology can be summarized as follows: 

• The formalism is based on the Polyakov world-line path integral, whose particular 
merit is that it is structured in terms of the spin factor, a quantity which is path- 
determined and which accounts geometrically for the spin of the propagating particle- 
entity. Within this framework, the quark-gluon dynamics are embodied in the expec- 
tation value of open/closed Wilson lines. Thus, quantities of interest for us, such as 



5 We have implicitly assumed the pre-confinement property, originally articulated in the first work 
of reference [25] (see also the second one), according to which the non-perturbative dynamics 
responsible for confinement screens color up to the infrared scale A which sets the lower limit for 
the perturbative regime. 
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amplitudes or cross sections can be solely described in terms of appropriately weighted 
integrals over Wilson contours. 

• We have based our considerations on the large scale Q 2 in order to apply the math- 
ematically well-founded fact that the local characteristics (like endpoints, cusps, etc.) 
of the contours involved in the path integral can be dissected out (cf. Eq. (52)). In 
this way we calculated the resummed expression for soft-gluon emission that gives rise 
to the Sudakov factor. 

• We reiterate that our primary goal has not been to reproduce known results, but 
to show how to obtain them on the cross-section level using the world-line Polyakov 
path integral. It is obvious that this type of approach can be used to describe DIS- 
type processes as well (see [26] from which the present investigation partly derives). 
Furthermore, since we are not obliged to use special type of paths and we can always 
stay within the Euclidean formulation - at least as long as we care about the leading 
behavior of quantities like the cross section, the present investigation may pave the way 
to extend this type of approach to the large transverse distance regime, where we shall 
meet power corrections signaling nonperturbative contributions [21,24,27-30]. Thus, 
we hope that this formalism can be extended to the calculation of the non-perturbative 
effects of the expectation value of Wilson loops by having recourse to the extensive 
existing literature [9,31] on the subject. 
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APPENDIX A 



Our task is to establish Eqs. (19) in the text. Performing the integration entering the 
right hand side of Eq. (18), one obtains 



r(«) 



1 

4^ 



x{„f(i,--1; 



2A (4-D)/2 r /D \ _l 1_ 

V2 / A-D D-3 
D . D-l 



2w 



-;l-w' 



D-l 1 + w 



4[2(l--)] 2 - D/2 ^(l,y-l; 
Setting D = 4, we obtain 

F(l,l;3/2;1-™ 2 ) = — - 

v y sinn 7 cosh 7 



7 



and 



F(l,l;3/2;- 



w 



7 



7T 



— I- 



sinh 7 sinh 7 



(Al) 



(A2) 



(A3) 



As the imaginary part in the above expression will cancel against its counterpart in the 
conjugate expression, it can be dropped as far as the cross-section is concerned. 

Denoting the expression inside the curly brackets on the rhs of Eq. (Al) by fniw) and 
setting 



f D {w) = h{w) + {A-D) 



fpjw) ~ U{w) 
A-D 



(A4) 



we realize that the second term on the rhs will lead to finite terms that depend solely on w 
and which will cancel against similar contributions of the same sort coming from the other 
terms entering Eq. (12). Putting everything together, one finally arrives at Eq. (19). 
To establish the result given by Eq. (22), we first note that Eq. (21) gives 



_}_(_ V_ 

An 2 [ 71 A 2 



2\ (4-0)/2 



2 J (A-D) 2 



2vv' 



\v\ 



(4-D)/2 



X 



„ (D , n D D 2vv'\ ( 2vv' 
F --1,2--; 3--; — ) + [ 1 + 



2-D/2 



2 2'~ 2' \v\ 2 \ \v\ 2 



2v ■ v' 



\v\ 



2-D/2 



(A5) 



Then, in the limit D — * 4 one easily retrieves Eq. (22). 
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